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Abstract. We propose a way to find the asymptotic distribution of zeros of 
orthogonal polynomials p n (x) satisfying a difference equation of the form 

B(x)p n (x + 5) - C(x, n)p„(x ) + D(x)p n (x - S) = 0. 

We calculate the asymptotic distribution of zeros and asymptotics of extreme 
zeros of the Meixner and Meixner-Pollaczek polynomials. The distribution of 
zeros of Meixner polynomials shows some delicate features. We indicate the 
relation of our approach to the WKB technique and to the approach based 
on the Nevai-Ullman distribution. 
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1 Introduction 


The hypergeometric orthogonal polynomials p n (x) that appear in the Askey 
scheme [jl] satisfy either differential equations or difference equations of the 
form: 

B(x)p n (x + 6) - C(x, n)p n (x) + D(x)p n (x - S) = 0, (1) 

where <5 is equal to 1 or i when the polynomials are orthogonal with a discrete 
or continuous weight a(x), respectively. 

As is well known (see, e.g., @), the zeros of p n {x) are real and simple. In 
the present paper, we will propose a procedure to calculate the asymptotic 
distribution of the zeros of p n (nz ) as n —> oo by applying a form of Bethe 
ansatz to equation (1).Q The difference between the cases 8 — i and 8 = 1 
will be clearly visible. 

Note that although we focus on the orthogonal polynomials belonging 
to the Askey scheme, our results have a more general application. In fact, 
we need only the two conditions: 1) polynomials p n (x) satisfy the difference 
equation (1) (or a certain type of (/-difference equation |J); 2) the density 
of zeros is piece-wise smooth on a one-dimensional manifold. For the poly¬ 
nomials in the Askey scheme, the latter manifold is just the real line. For 
another example, where the manifold is the unit circle, see |3j]. 

In this work, we are interested in zeros of p n (x) only to the main order 
in n. However, the method we apply should remain useful for extraction of 
the further terms in the asymptotic series. 

Other approaches to the calculation of the asymptotic distribution of 
zeros of orthogonal polynomials are based on the knowledge of the recurrence 
relation, or the weight function da(x)/dx, or the asymptotic behaviour of 
polynomials themselves (see [fij for a review). 

Our approach is closely related to the semiclassical (WKB) analysis of 
differential (difference) equations. WKB techniques easily yield the more 
“rough” features of the distribution of zeros such as their density. Through¬ 
out the paper, we indicate how some of our results can be obtained by the 
WKB method. 

We will consider in detail the cases of the Meixner polynomials (orthog¬ 
onal on the set {0,1,2,...}, (5 = 1) and the Meixner-Pollaczek polynomials 
(orthogonal on the set (— 00 , 00 ), 5 = i). The asymptotics of zeros of the 

1 We will consider infinite systems of orthogonal polynomials. 
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Meixner-Pollaczek polynomials has been investigated in |||7| by other meth¬ 
ods. Therefore, the reader can compare several approaches. 

The zeros of p n (x) are intimately related to the points of increase of 
the weight (see, e.g., ||). If the moment problem corresponding to the 
polynomials p n (x) is determinate, the points of increase of the weight coincide 
with the closure of the set of zeros of p n (x), n —► oo. For the Meixner and the 
Meixner-Pollaczek polynomials the latter fact implies that the density p(z) 
of the zeros of p n (nz) at the point z = 0 is 1 and oo, respectively. Naturally, 
the expressions for p{z) we obtain have this property. It is interesting that, 
moreover, for the Meixner polynomials p(z) = 1 in the finite vicinity of 0 = 0. 
In this vicinity the distribution of zeros is characterized by another function 
(simply related to p(z)) that gives the rate with which the distance between 
zeros approaches the value 1 as z —> 0. This rate is exponential in n. 

There is another interpretation of our results on the Meixner and Meixner- 
Pollaczek polynomials. Consider a symmetric tridiagonal (Jacobi) matrix 
with real matrix elements (= Jki, i, k — 0,1, 2,..., Jik = 0 if \i — k\ > 1) 
such that the limits lim ?woo J nn /<^(n) = a, lim^oo J n n+ i/(p(n) = 6/2 are 
finite and not simultaneously zero and <p(n) is a nondecreasing function. All 
such matrices for tp(n) = n are separated in two classes. The first class is 
comprised of matrices related to the Meixner polynomials and the second, to 
Meixner-Pollaczek. The asymptotics of zeros of polynomials is the asymp¬ 
totics of eigenvalues of the corresponding truncated Jacobi matrices. Thus 
in this paper we have calculated, in particular, the asymptotic distribution 
of eigenvalues of the above described Jacobi matrices for <p{ri) = n. On the 
other hand, the density of zeros in the case of a = 0 is known for any <p(n) 
(with some limitations on the form of ip(n)). It is called Nevai-Ullman dis¬ 
tribution 1JH- It can be easily generalized to the case of any real a, b. In 
particular for <p(n) = n, this gives an alternative derivation of the density of 
zeros of the Meixner and Meixner-Pollaczek polynomials. However, this does 
not reproduce the delicate features of the zero distribution of the Meixner 
polynomials in the vicinity of z = 0. 

The plan of this work is the following. I 11 Section 1 we introduce the Bethe 
ansatz technique in the simplest case of polynomials satisfying differential 
equations. We obtain the asymptotic density of zeros and compare the results 
with those of a more straightforward approach. 

In Section 2 we consider the difference equation (1) with S = i. It is satis¬ 
fied by polynomials orthogonal with a continuous weight: Meixner-Pollaczek, 
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continuous Hahn, continuous dual Hahn, Wilson. After getting the general 
solution, we focus on Meixner-Pollaczek as the other polynomials have, to the 
main order in n, the density of zeros simpler than that of Meixner-Pollaczek. 

In Section 3 we consider the more interesting case of equation (1) with 
5 = 1. Such equation is satisfied by the Meixner and Charlier polynomials. 
We will consider only them, as other polynomials with discrete weights in 
the Askey scheme are finite systems. 

In Section 4 we consider the asymptotic eigenvalue problem for symmetric 
tridiagonal matrices. 


2 Differential equation 


We shall start with the case of a differential equation as it is the simplest 
one, and it already shows some features of the method. 

Consider a system of polynomials satisfying the following differential 
equation:^ 

a(x)y'n + K x )y' n + c(z, n)y n = 0. (2) 

Here y n {x) is a polynomial of exactly degree n, functions a(x), b(x), c(x,n), 
where only the last one depends on n, are smooth with possible singular 
points that do not coincide with any zeros of p n (x). We shall assume that 
there is such a constant y > 0 that c(x,n) ~ x )n 2/i as n —> oo. Such 
equations are satisfied, e.g., by the Jacobi, Laguerre, Hermite polynomials. 
These polynomials are orthogonal with an absolutely continuous weight on 
an interval A (unbounded in the case of Laguerre and Hermite) and the zeros 
of y n (x) become everywhere dense in A as n —> oo. 

Let x\ < X 2 < • • • < x n be the zeros of y n {x). Up to a constant factor, 
we can write y n (x) in the form y n (x) = rifc=i — x k)- Substituting this 
expression in (S[) and dividing the equation by y n ( x ), we get 


a( x ) E E 

k=1 i=1 


1 

(x - x k )(x 


x i) 


n i 

+ b(x) E-b c(x, n ) = 0. 

ife l x ~ x k 


( 3 ) 


Obviously, the singularities of the function at the left-hand side must vanish. 
Equating the residues in the simple poles x k , k — 1,..., n, to zero, we obtain 

2 Here we restrict our attention to second-order differential equations, but it is easy to 
generalize our approach to higher-order equations. 
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the following system: 


2 a(x k ) E 


i=i %k %i 
i^k 


+ b[xk) = 0, k = 1,... ,7i. 


( 4 ) 


Equations (fl) can also be obtained by dividing (J2j) by y' n (x) and putting 
x = Xk- The quantities a p k = Y^= i 1 /(xk — Xi) p , p = 1, 2,... and related ones 

i^k 

were calculated in |[T2| . [TTj| . [T4|| . Analogous algebraic relations for solutions of 
some g-difference equations are given in [15]. Because of the analogy to the 
theory of integrable models, the equations of type (f|) are called Bethe ansatz 
equations. For our purposes, it will be sufficient to know only the asymptotics 
of aik and a 2 k as n —> oo. 

Let us now differentiate (H) with respect to x and apply the same polynomial- 
ansatz technique to the resulting equation. Then instead of ( 0 ) we get 


3 a(x k ) < J2 


E 


-E 


+ 


i =i Xk Xi j =1 Xk Xj i=1 (Xk Xi ) 

, i^k i^k 

n ^ 

2 (a'(x k ) + b(x k )) E-1- c ( x k, n) + b'(x k ) = 0, k = l,...,n. 

~i Xk - Xi 

i^k 


( 5 ) 


Now divide (||) by ri 2p and take the limit as n —> oo on condition that the 
point x = lirn^oo Xfc( n ), x G A, is fixed. Then, using the fact that, as follows 
from (fl), 


1 


lim — 

n .— VfYl ry-iCT t ^ 


n—>00 nri® - -- qp - rp . 

11 i= 1 

i^k 


= 0 


( 6 ) 


for any a > 0, we obtain from 

1 


i. \ ^ 1 Coo(x) 

lim —— > --— = ——— 

n^oo n ^ i^ Xk - Xi y 3 a(x) 

i^k 


( 7 ) 


Let us assume that the zeros are asymptotically spaced as l/n M and their 
asymptotic density is a piece-wise smooth function. Hence for \i — k\ < M 
we have 

Xi - x k = (i - fy'fk/n* 1 (8) 
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up to at most 0((M/n M ) 2 ) as M,n —> oo in such a way that M/n M —> 0. 
For brevity, we shall denote the limit M,n —> cx), M/n M > 0 just by lim. 
In (j8j) lim 7 fc(„) = 7 ( 0 ;) > 0. The asymptotic density of zeros is, obviously, 
p(x ) = lim l/ 7 fc ( n ). Using (||) we rewrite the l.h.s. of (0) in the form: 


lim 


1 


E 


n 2fi ^ _ x .y 

i^k 


00 1 7 r 2 

= 2p(x) 2 E- = -^-p(ar) 

i=i ? J 


( 9 ) 


So putting (^) and (|9|) together, we have 


p(x) 


1 Cqq{x) 

7 a(x) 


( 10 ) 


on the set A. It is not necessary to justify ([h]) further as this result is, of 
course, well known. It is related to the semiclassical approximation in the 
theory of Schrodinger equation (see, e.g., [|I7]]). In a few words, the usual way 
to obtain (|T0|) from (^j) is the following. Consider a small cr-neighbourhood of 
a point x. If a(x)c 00 (x) > 0 then, starting with a sufficiently large number n, 
the characteristic values of the equation (||) (considered as an equation with 
constant coefficients: a(x), b(x), c(x,n ) being sufficiently smooth), u± = 
{— b(x) ± \Jb(x) 2 — 4a(x)c(x, n)}/(2a(x)), have the imaginary part growing 
as n Hence, it is easy to conclude that the asymptotic distance between 

consecutive zeros is Xk+i — Xk r\_; 7fc = 7rTE^+o(l). 


For the polynomials orthogonal on an unbounded interval, another type 
of the density of zeros can be introduced. Take, for example, the Hermite 
polynomials. Their zeros are symmetric with respect to the origin, and the 
largest (smallest) zero x n ~ \/2n (aq ~ — y/2n). Equation (£|) in this case 
has the form: 


n y 

E-= x k , k = l,...,n. (11) 

Xk ~ x t 

i^k 


Changing the variable aq = z i \f 7 2n and then taking the limit n —> 00 , we get 


0 

ly f 1 P(w)dM 
2 J- 1 z — LU 


Z e (-1,1), 


( 12 ) 
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where the density of zeros p{oj) = lim^oo l/(z i+ i — Zi)n, u; = lim^oo Zi( n ). 
Solving (p~2|) and using the normalization condition f\p{uj)duj = 1, we get 


p(z ) = —Vl — z 2 . 

IT 


(13) 


This way to derive (|T3|) was proposed in 
to find p(z). Expression ([ 


If 


Of course, there are other ways 
is valid for a contracted to [—1,1] density of ze¬ 


ros of a large class of orthogonal polynomials (to which Hermite polynomials 
belong), ft is a particular case of the Nevai-Ullman distribution (see Section 
5) for p(n) = sjn. 

Note that (|10j) for Hermite polynomials gives p(x) = \/2/tt. 

In the next Section we shall see how the technique we used to obtain ([10) 
generalizes to the difference equations. 


3 Difference equation. Continuous weight. 


Consider polynomials in the Askey scheme orhtogonal with a continuous 
weight on an unbounded interval and satisfying (1). In this case 5 = i. 
Take the n’th polynomial and make the change of variable x = zn. Then 
Pn{x) = nLi(® - x k ) = n n y(z ), where y{z) = TILi ( z ~ z k )■ before, 
we assume X\ < x 2 < ■ ■ ■ < x n . If C(x,n ) ~ c(^)n M as n —> cxo, then set 
b(z ) = lim B(z7i)/ ji 11 , d(z) = lim^oo D(zn)/n fl . Thus, to the main order 
in 7i, (1) can be written in the form: 

b(z)y(z + i/n) — c(z)y(z) + d(z)y(z — i/n ) = 0. (14) 


Substituting here y(z) = YYk=i ( z ~ z k ) and evaluating the equation at the 
zeros z m , m — 1,2, , n , we get: 


n 


k =1 


z m z k + i/n 

Z m - z k - i/n 


diyZm) 

b(Zm) ’ 


m — 1,2,..., 7i. 


(15) 


Now assume that the distance between zeros Xk is of order 1 (which is 
natural to expect from (1)), and the asymptotic density p(z m ) = l/{x m — 
x m -i) of zeros becomes a piece-wise smooth function of z — lirr^^oc z m ( n ) as 
n —> oo. Hence for \p\ < M 


Z m —p 


P 


p(z)n 


(16) 
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up to at most 0(M 2 /n 2 ) as M, n, n/M —> cx). As before, we denote this limit 
by “lim”. Applying “lim” to (|T5|) we have 


lim 


/ M -h + i + i 

-n ^_ 8 ^_, n 

\ P =1 p(z) L p(z) k 


1 + 
1 - 


i / n \ 
z m Zp | 

i / n I 

Zm — Zk / 


d(z ) 

WY 


(17) 


where the prime indicates that the product is taken over k outside the range 
m — M ,..., m + M. The product over p gives 1 in the limit. Since i/n(z m — 
Zk) in the second product is small (of order 1/M or less), we can perform 
the small-parameter expansion in it. Replacing the second product by the 
exponent of its logarithm we get limexp[(2i/n) Xb- l/(zm ~ z k)\ = d(z)/b(z) 
which, in turn, gives 


exp 


— 2i V.p. 


r°° p{uj)duj N 


1-00 iv — z 


d(z) 

b(z)' 


(18) 


Note that p(z) also satishes the normalization condition 



1. 


(19) 


Instead of trying to solve (|1^) directly, we will do the following. Replace 
z by z -\- i/n in (|T4|) and consider this new equation (cf. Section 2). Then we 
obtain instead of (pA|): 


n 


z m z k T 2?/ 11 c(xO m ) 


k=l A " m Ac A ^(^m) 

Proceeding as before, we get for the l.h.s. of 


m = 1 , 2 , 


,n. 


( 20 ) 


z m Zk T 2% j11 


[M/2] 


lim TT —-— -= 2 lim TT { i + 

k=i An - z k + i/n { 

p(iv)dui\ 


' 2 P {z) ' 
2s + 1 


exp 


Zk, 


2 cosh(7rp(z)) exp ( — i V.p. J 


— OO CO — Z 


( 21 ) 


Thus, in the limit equation (|20|) becomes 


r°° p(u)d(v\ c(z) 
) b(z) 


2 cosh(7rp(z)) exp — iV.p. 


— OO UJ — Z 


( 22 ) 



























Putting (f[|) and (|22|) together, we obtain 


p(z) = — arccosh 

7T 


c{z) 


2 \Jb(z)d(z) 


(23) 


Thus p(^) can be either (|23|) or zero. If we suppose that p(z) is defined by 
(|23|) for all z where \c(z)/(2^jb(z)d(z))\ > 1 (hence p(z) = 0 otherwise) then 
the asymptotic values of the largest and smallest zeros are obtained from the 
equations c(z) / (2^b(z)d(z)) = ±1. 

In the concrete cases we can verify if thus obtained p(z) is indeed the 
asymptotic density of zeros by substituting it in (|l8|), ( |l9|) , (^) and integrat¬ 
ing. These concrete cases include the Wilson, continuous Hahn, continuous 
dual Hahn, and Meixner-Pollaczek polynomials |l|. 

The most interesting example is provided by the Meixner-Pollaczek poly¬ 
nomials. They are orthogonal on the line (—oo, oo) with the weight function 
da(x)/dx = exp[(2 <f> — 7r)x]|r(A + ix) | 2 , where A > 0, 0 < <f> < it. (Note that 
dafx) can be obtained H from the difference equation). The coefficients in 
equation (1) for these polynomials are the following: B(x) = — ix), 

C(x, n) = 2 i[(n + A) sin d — x cos0], D(x) = + ix), and 5 — i. 

Using the just obtained general results, we can immediately formulate a 
theorem. 

Theorem 1. For the Meixner-Pollaczek polynomials p n (zn) the asymptotic 
density of zeros is 


p(z) = —arccosh 

7T 


and p(z) = 0 otherwise. 


Sill i 


— COS ( 


if z G [— cot —, tan-^], 


We note that the values ntan | and —ncot | are the asymptotics of the 
largest and smallest zeros of the Meixner-Pollaczek polynomials p n (x), re¬ 
spectively. There can be no isolated zeros z m outside the interval [— cot |, tan |] 


as the Bethe ansatz equations (p~o| . p0| ) would not hold for such zeros. We can 
also use here another argument. Consider the Gerschgorin circles of the Ja¬ 
cobi matrix associated with p n (x) (see Section 5). We then find that all zeros 


must be inside the interval [— cot |, tan |]. 
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On the assumption of smoothness of p(z), it is easy to prove the theorem 


by verifying (18), (19), and (22). In particular, we have 


V.p. 


f tan I p(co)duj [0 — 7T, 0 < z < tan■ 


/— cot ^ UJ — Z 


— cot | < z < 0 


(24) 


As such, our results for the Meixner-Pollaczek polynomials are not new. 
In the implicit form they are contained in 

In the next section we will use Bethe-ansatz techniques to indicate some 
interesting features of the distribution of zeros of Meixner polynomials in 
the region where the asymptotic density of zeros is constant. However, the 
formulas for the density itself ((|25|) and (|33|)) are easier to obtain using the 
WKB approach which we now briefly describe. Suppose that y(z) has the 
asymptotic form y(z) ~ exp(nS'o(^) + S'i(^) + S , 2 (^)/n + - • •), where Sj(z), j = 
0,1,... are continuous, and substitute this expression into (O). Decompose 
Sj(z ± i/n) into Taylor series and collect the terms of the same order in n 
in the equation. Solution of the equation to the main order in n gives S 0 (z) 
which, in turn, yields the density of zeros. We see that this approach gives 
not only the density of zeros but also the asymptotic expression for y(z). 


4 Difference equation. Discrete weight. 

The Meixner and Charlier polynomials orhtogonal on the set {0,1,2,...} 
satisfy equation (1) with 8 — 1. Using the same notation as in Section 3, we 
reduce (1) to the form: 

b(z)y(z + 1/n) — c(z)y(z) + d(z)y(z — 1/n) = 0. (25) 


Evaluating it at the zeros z m , m — 1, 2,..., n, we get: 

tt z m T 1/n d(^ m ) 

11- TT = ~T7 —V’ m = 1 

k=! z m- Zk- 1/n b(z m ) 


(26) 


For the distance between nearest zeros we still have the expression flTBj ). 
However, as we can notice from (^6j), there appears a new important feature. 
We mentioned in the introduction that in the present case p(z) —>• 1 as z —> 0. 
Hence we can expect that z m — z m -\ — l/n = o{l/n) and z m — z m+ 1 + 1/n = 
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o(l/n) for z m close to zero. Therefore, to the main order in n, the ratio of 
these quantities is an unknown if p(z m ) = 1 and equal to —1 if p(z m ) ^ 1.0 
We denote A m -i = z m — z m -i — 1/n. Now we can proceed as in the previous 
section. 


.. tt z m z k T 1/n 

inn -—— = inn 


Zm - Zm+ 1 + 1/n Zm ~ Z m -1 + 1/n 


k =1 


Zm- z k - 1/n 


Zm 

1/n 


M _E_ + 1 _z£_ i 1 1 I 

n P(z) P( z ) TT' Zm~Z k 

P 1 -p i 11 , 


JP _I JZP. _ 1 

P = 2 p(z) p(z) k 


lim 


A r 


, A m _i 


exp 


2 , I 

n k z ” 


1/n 
Zm ■Z/c 

1 


Zm+l - 1/n Zm - Zrn-1 “ 1/n 


d(z) 


x 




Thus, in the limit fl26|) becomes 
X(z) exp -2 V.p. 


Hz)' 


r°° p(co)duj\ d{z) 


l—oo (jJ — Z 


Hz )’ 


where x(z) = lirn(A m /A m _ 1 ), z = lim^oo z m {n)■ 

Remember that p(z) satisfies the normalization condition 


/ OO 

p{z)dz = 1. 

-OO 


(27) 


(28) 


(29) 


As in Section 3, change the variable z —> z+l/n in d25|) . Evaluating thus 
obtained equation at the zeros z m , we have: 


tt z m zk T 2/n _ c(z m ) 

til z m - Z k + 1/n b(z m ) ’ 


(30) 


In the limit ([30|) becomes 



(31) 


3 Strictly speaking, we also had to consider the cases when p(z m ) = s, where s is an 
integer greater that 1. In this case z m — z m ± s ± 1/n = o(l/n). However, we shall see that 
p(z) < 1 for all z at least for the Meixner and Charlier polynomials. 
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Since lirn(A m+ i/A m ) = x( z )i we may rewrite ([H]) in the form: 

r°° p(co)duj\ c(z ) 


(1 + x( z )) cos(vrp(^)) exp | -V.p. J 
From (|28|) and (E 


-oo U) — Z 


b(z)' 


( 32 ) 


we obtain a simple relation between p(z) and x( z ) : 
c(z) x( z )b(z) 


(1 + x(z)) cos(vr p(z)) = 


b(z ) \ d(z) 


(33) 


We know that x( z ) = 1 whenever p(z) is nonconstant. Hence, if p(z) is 
nonconstant, it satisfies the equation: 


cos(7rp(2:)) = 


c(z) b(z) 


(34) 


2 b(z) \ d(z) 

Obviously, this equation can hold only in the region where the absolute value 
of the r.li.s. is no more than 1. Having guessed p(z) and x( z ) from (|33|) , 
we can check if we are right by verifying (|28|), (|32D, and (^). Thus we 
get the following theorem for the Meixner polynomials p n (x ) for which [jlj 
B(x) = ( x + /3)c , C(x) = x+ ( x + /3)c+ (c— l)n, D(x) = x, (3 > 0, 0 < c < 1. 
Theorem 2. For the Meixner polynomials p n (zn) the asymptotic density of 
zeros is 

1, z G [0, a -1 ] 

p(z) = \ ^arccos f(z), z e [cc -1 , a] , 

0, z [0, a] 


where 0 < arccos z < 7T, and 

1 + \Jc 


a = 


f( z ) = 


(c+ l)z + c — 1 
2 z/cz 


The function 


In 


vio = / 2arcosh \ f( z )l Z e[0,a x ] 
X{ * } \0, ze[a-\a]' 


Remarks. Since in the limit A m _i = exp[— (k/n)nlnx(z)]A m -i + k, and 
x{z) > 1 for z G [0,cc _1 ], the function lny(^) describes the exponential rate 
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with which zeros Xk approach integer values as z decreases from the point 
ct _1 . 

The function p(z) is continuous on the set (0, oo). Its derivative, however, 
is discontinuous at the points a -1 and a. 

The value is the asymptotics of the largest zero x max of the Meixner 

polynomials p n (x). Note that the inequality x max < + °{ n ) is required 

by the position of the Gerschgorin circles of the associated Jacobi matrix. 
For further and more precise inequalities for zeros of p n (x) see [§]. 

By the substitution 0 = In ^/c, we can make similarities with the Theorem 
1 more transparent. Indeed, f(z) = (sinh0)/z + cosh0, a = — coth(0/2). 
One might have expected such similarities, because the Meixner-Pollaczek 
polynomials are obtained from the Meixner polynomials by continuing the 
latter in the variable x and parameter c into the complex plane ||. 

The asymptotic equation (f25[) for the Charlicr polynomials can be re¬ 
garded as a particular case of the one for the Meixner polynomials when 
c —* 0 . 

In the next section we will consider another approach (initiated in |J) 
to calculation of the asymptotic density of zeros of orthogonal polynomials. 
It is based on the study of the recurrence relation rather that the difference 
equation. This approach allows us to give an alternative derivation of p(z) 
for Meixner and Meixner-Pollaczek polynomials (Theorem 1 and the first 
part of Theorem 2). 


5 Eigenvalue density of Jacobi matrices 

Consider a symmetric Jacobi matrix J with real matrix elements (= Jki , 
i,k — 0,1,2,..., J ik = 0 if \i — k\ > 1) such that J nn ~ a</?(n), J nn + i ~ 
hp{n )/2 as n —> oo, where a and b are not simultaneously zero. The function 
<p(x) : R + -> R + is nondecreasing and such that lim^oo <p(n + x)/ip(n) = 1 
for any x G R; also let lim )WOO ip(nt)/(p(n) = ijj(t) exist for t € [0,1], and the 
function g(t) = dt/d'ip be continuous. We will call the class of such matrices 

Let J{n) be the (n + 1) x (n + 1) truncated J\ J(n) = \\Jik\\ik=o- ^ 
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follows from [p|,[TTj| that for k — 0,1,.. .| 

k 


p k = lim — Tr 

n—>oo fi 


J{n) 

<p(n) 


1 r 1 , r a +\ b \ 

— / i/ k {t)dt 
7T JO 


x k dx 


°-l b l Jb 2 — (x — a)' 2 


( 35 ) 


Suppose there exists a solution p(z) to the moment problem p k = z k p(z)dz, 
k = 0,1,.... Since a and b are finite, the function p(z) = p(z,a,b ) has a 
bounded support which implies that the moment problem is determinate, 
and hence, p k , k — 0,1,... uniquely define p(z, a, b). The function p(z, a, b) 
is the density of eigenvalues of J(n)/<p(n) in the limit n —»■ 00 . So up to a 
uniform scaling of the spectrum, p(z, a, b) depends only on a/b. 

On the other hand, recall that the orthogonal polynomials satisfy the 
recurrence relation [2[: 


XPn( x ) — a n p n +l + (3 n Pn + InPn-li 71 = 0 , 1 ,..., 7oP-l = 0 . ( 36 ) 

After the transformation q 0 = Po, q n {x) = \JPn(x) , n = 1 , 2 ,..., it 
takes the form 

•^'9re(*7') Jnn+lQn+1 T JnnQn T Jnn— lQn— 1) 77. 0, 1, . . . 

Jnn Jnn+1 'A+l n \/Ctnjn+l> To — lQ— 1 0, 

where J = | Afc| |ffc=o is a symmetric Jacobi matrix. As is well known (e.g., 
[0), the zeros of pn + i(x) are eigenvalues of J(N) = || Jik\\ik=o- Thus the 
asymptotics of the zeros of p n (x) gives the asymptotics of the eigenvalues 
of J(n). Notice, that if J E A^ r p and p(z,a,b ) is defined as above, we get 
p(z, a, b) = p(z, a, —b) by the substitution q n (x) = (—1 ) n q n (x). 

Now let us consider the symmetric Jacobi matrices associated with the 
Meixner and Meixner-Pollaczek polynomials. For Meixner-Pollaczek ||T[ we 
get 


2 sin < 


\ yit ~r J-y i 

Jn n+1 = Z T T 

2 sm 0 

Thus, J E A n and a/b = — cos (j), (j) E (0, n). So —1 < a/b < 1. 

4 The values n k for 6 = 0 are obtained from (|35|) using continuity in b. 
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J 


For the Meixner polynomials (0 = In y/c) 

\J i n + l)( n + fi) c n 


n n+1 


C — 1 


2 sinli 0 ’ 


Jnn 


n + (n + 0 )c 

1 — c 


—n cotli 0. 
(39) 

Thus, J £ A n and a /6 = — cosh0, 0 G (—oo,0). So —oo < a /6 < —1. 

Now it is easy to verify that, up to a sign, an arbitrary matrix J from A n 
for which |a| ^ \b\ has the form r Jo, where r = J\a 2 — b 2 | and Jo is a matrix 
associated with either Meixner-Pollaczek or Meixner polynomials defined up 
to a sign. Obviously, p(z,a,b ) = po(z/r,a/r,b/r)/r , where p and po are the 
eigenvalue densities of J(n)/n and Jo(n)/n as n oc, respectively. 

Now we can reformulate Theorem 1 and the first part of Theorem 2. 
Theorem 3 . Let J G A n , J nn ~ an, J nn +i ~ bn/2 as n —> 00. T 6 en i/ie 
eigenvalue density of J{n)/n in the limit n —► 00 zs following: 

1) 0 < a /6 < 1 , 

^ G [a — 6 , a + 6 ] 

2 ^ [a — 6 , a + 6 ] 



) —arccosh 

/ nr 

rl 4. °l 

zb ' b 

lo, 



P{z,a,b) = 


i/r, 

ijarccos (-g + f) , 

0, 


z G [0, a — 6 ] 

2: G [a — 6 , a + 6 ] , 
2 0 [ 0 , a + 6 ] 


where 0 < arccos((c) <7 t, r — J\a 2 — 6 2 |. 

In the case a/6 > 1, more precise information about the zeros in the 
region where p(z ) = 1/r is given by the function x{ z ) (Theorem 2). 

It was sufficent to consider only the case a > 0, 6 > 0 (substitution 
6 —> — 6 does not affect the spectrum, and a —> —a inverts it with respect to 
zero). I 11 what follows, we always assume a and 6 nonnegative. 

For a = 0 our result coincides with a particular case of the Nevai-Ullman 
distribution. This distribution is the following fll0| : 

1 f 1 gfjj)du} 


p(2, 0,6) = - f 

IT J\z 


7T J\z/b\ y/b 2 UJ 2 - 2 2 

For <p(n) = n 1 , 7 > 0, we have g(u>) = ui~ l+l ^ /^^\ 
5 The results for 7 = 0 are obtained using continuity. 


2 G [—6,6]. 


(40) 
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It is not difficult to generalize to arbitrary a following the approach of 
0 .0- It is, however, interesting because of qualitatively new features which 
appear when a > b. Note that this generalization appears for <p(ri) = n 7 in 


11], where a different approach to its derivation is used. 


Theorem 4. Let J G J nn ~ cup(n), J nn +i ~ btp(n)/2 as n —> oo. 

Then the eigenvalue density of J(n)/(p(n) in the limit n —► oo is the following: 
1) 0 < a/b < 1, 

pO) = \ fz/(a+bsi gnz ) ^ b2 °^_ auj) 2 l f z G [a~b,a+b], and p{z) = 0 otherwise; 


2) a/b > 1, 

P{z) = 




+<U/<* + »> yp5+=?. *e[0,a-6] 

1 £/(«+*) ze{a-b,a + b]- 

0, z [0, a + b] 


Note that for <p(n) = n 7 , 7 > 0 we can write case 2 in a simpler form: 

,-1+1/7 


,, /(«-»)(*)' 

P( z ) = \ h(z), 

0 , 


where 


1 r 1 

h(z) = — / 

7ry Jz/ 


, z G [0, a — b] 
z G [a — 6, a + b] > 
z [0,a + 6] 

to - i + i/ 7 (ia; 


7T7 J+(a+ft) +2^2 _ (~ _ aw )2 


Proof is analogous to that given by Ullman in the case a = 0. Rewrite 
fl55|) in the form: 


Pk 


1 r 1 


7r jo 




r a+fe ( ljjx) k dx 


a ~ b yb 2 ~ (x — af 


(41) 


Consider case 1. Changing the variables z = ipx, uj — if in the double 
integral, we have: 


Pk 


i rv* r 


g(cu)dw 


7T J a—b 


Jz/{a+b sign z) ^J}) 2 UJ 2 — (z — OC<j) 2 

which proves case 1. Similarly, we prove case 2. 


( 42 ) 
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Theorem 1 and the first part of Theorem 2 (that is Theorem 3) are corol¬ 
laries of Theorem 4 when <p{n) = n. 

Finally, note an important and, in its generality, very difficult problem: 
what weights of orthogonal polynomials lead to the matrices in A 

recent review of related results is given in [[J]. In this sense, the case |a| > |6| 
is studied less than the case |a| < |6|. 
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